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Abstract 

An effective method for ultrafast magnetization reversal of nanoclusters is sug- 
gested. The method is based on coupling a nanocluster to a resonant electric circuit. 
This coupling causes the appearance of a magnetic feedback field acting on the clus- 
ter, which drastically shortens the magnetization reversal time. The influence of the 
resonator properties, nanocluster parameters, and external fields on the magnetiza- 
tion dynamics and reversal time is analyzed. The magnetization reversal time can be 
made many orders shorter than the natural relaxation time. The reversal is studied for 
both the cases of a single nanocluster as well as for the system of many nanoclusters 
interacting through dipole forces. 
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1 Introduction 



The effect of magnetization reversal in nanomaterials is of considerable importance for vari- 
ous magucto-electronic devices, magnetic recording and storage, and other information pro- 
cessing techniques. The standard way of recording an information bit is to reverse the 
magnetization by applying a magnetic field antiparallel to the magnetization. From another 
side, the nanoparticlc magnetic moment has to be sufficiently stable, which can be achieved 
by the use of materials with high magnetic anisotropy. But the latter complicates the process 
of magnetization reversal. In order to resolve the contradiction between these two require- 
ments, different methods of magnetization reversal have been suggested, as can be inferred 
from the review articles [1,2]. 

Magnetization reversals caused by thermal fiuctuations and phonon-assisted quantum 
tunneling are rather slow processes at low temperatures, below the blocking temperature, 
where nanoclusters exhibit stable magnetization [1,2,3-10]. To make the reversal faster, sev- 
eral methods have been suggested. Thus, one can employ transverse magnetic constant fields 
[11] or short pulses [12-17], transverse microwave alternating fields at magnetic resonance 
frequency [18-24], and optical laser pulses [25]. 

In the present paper, we suggest another method for achieving an ultrafast magnetization 
reversal of nanoclusters. The method is based on coupling the considered nanocluster with 
a resonator by placing the cluster into the magnetic coil of an electric circuit. Then the 
motion of the cluster magnetization produces a magnetic feedback field acting on the cluster. 
This feedback mechanism essentially accelerates the magnetization reversal. Actually, the 
effect of the accelerated thermalization of nuclear magnets was proposed by Purcell [26] and 
considered, using the classical Bloch equations, by Blombergen and Pound [27]. Here we 
study the magnetization reversal by employing quantum microscopic Hamiltonians typical 
of strongly anisotropic nanoclusters possessing large spins. 

We study magnetic clusters with the effective sizes shorter than the exchange length of 
atoms composing them, when the magnetic cluster is in a single-domain state and its magne- 
tization can be represented by a large total spin. Such clusters are necessarily of nanosizes, 
which explains their name of nanoclusters. To avoid complications, due to distributions 
of particle sizes and shapes, we consider the magnetization dynamics of similar nanoclus- 
ters. There are two essentially different cases. One is the magnetization reversal of a single 
nanocluster. And the other is the magnetization dynamics in an ensemble of nanoclusters 
interacting through dipolar forces. We study both these limiting cases. 

2 Spin Dynamics of Magnetic Nanoclusters 

Let us, first, consider a single nanocluster. The typical Hamiltonian of a nanocluster, with 
the total spin S, can be written in the form 

H = -/ioB -S-D {S'f + D2 {S^f , (1) 

where //q = —^75 is the cluster magnetic moment, 75 ~ 2//^/^ is the gyromagnetic ratio, 
Bohr magneton, B is the total magnetic field acting on the cluster, D and D2 are anisotropy 
constants. This Hamiltonian reminds the classical Stoner-Wohlfarth model [28,29]. However, 
we start with the microscopic quantum Hamiltonian (1), where the spin vector is treated as 
an operator. This will allow us to explicitly define all system parameters and to take into 
account quantum effects that can be important for the dynamics of nanocluster assemblies. 



One usually represents the cluster energy in a reduced form [1,2] with the anisotropy 
parameters related to D and D2 as 



with Vi being the single-cluster volume and 5", the cluster spin value. The second-order 
magnetic anisotropy is caused by magnetocrystalline anisotropy, shape anisotropy, and sur- 
face anisotropy [1,2]. Sometimes, one includes the fourth-order and sixth-order anisotropy. 
However such higher-order anisotropy terms arc usually much smaller than the second-order 
terms, so their inclusion does not essentially change the overall picture. 
The total magnetic field is the sum 

B = Boe^ + He^ + BiBy (3) 

of an external constant magnetic field, generally having the longitudinal, Bq, and transverse, 
Bi, components, and of the resonator feedback field H directed along the axis of the coil, 
where the cluster is inserted to. The field Bq is directed opposite to the initial cluster magne- 
tization. The magnetic field, created by the coil, is described by the Kirchhoff equation. The 
latter, as is known, defines electric current generated in the coil by varying magnetization. 
In turn, the current produces magnetic field along the coil axis. The resulting equation for 
the generated magnetic field H — H{t) can be written [30,31] in the form 

^ + 2lH + co' £ Hit') dt' = ^ , (4) 

in which cu is the circuit natural frequency, 7 is the circuit attenuation, rj = Vi/Vres is the 
filling factor, with Vres being the volume of the resonant coil, and 

^ ^ (S^) (5) 

is the average magnetization density, corresponding to the mean spin in the direction of the 
coil axis. Thus, moving spins produce magnetic field that acts back on spins, accelerating 
their motion. 

To write down the equations of motion for the spin operators, it is convenient to introduce 
several notations. We define the Zeeman frequency 



the transverse frequency 



a;o = - y -Bo = 75S0 , (6) 



00, = -f^B, = ^sB, , (7) 



and the effective transverse field acting on the spin, 

f = -^lioH + uj,. (8) 
Then the Heisenberg equations of motion for the spin operators 5'='= = -S"^ ± take the form 

^ = -iuoS- + fS' + i^ (S'S' + S'S-) + i ^ (S'^S' + S'S'^) , (9) 
dt n n 



plus its Hermitian conjugate. And the equation for the longitudinal spin becomes 

= - [f*S- + S+f) + {S'^S^ + 5^5") . (10) 

UjL ft 

The measurable quantities are the average spin components, such as the reduced trans- 
verse component 

u^\{S-) (11) 
and the reduced longitudinal magnetization 

s ^ \ {S') . (12) 



Averaging the equations of motion (9) and (10), we use the decoupling 



(13) 



which preserves the correct behavior for all spins. Thus, for spin one-half it gives exactly 
zero, because of the anticommutativity of different Pauli matrices, and it becomes asymp- 
totically exact for large spins [32-34]. We also take into account that nanoclusters possess a 
longitudinal relaxation rate 71 that is due to phonon- assisted tunneling and to spin-phonon 
interactions of the nanoclustcr with its surrounding. 

Let us introduce the longitudinal anisotropy frequency 

un = (2S-1)^, (14) 

the transverse anisotropy frequency 

UJ2 = {2S-1)^, (15) 

the effective anisotropy frequency 

UJA = UJD + ^i^2 , (16) 

and the effective rotation frequency 

= UJq — ujas . (17) 

Also, let us define the effective field 

F = f + ^ UJ2U* = - ^ i^oH + ^^1 + ^ i^2U* . (18) 

Then averaging Eqs. (9) and (10) yields the equations for the transverse spin component, 

^ = -iflu + Fs , (19) 
at 

plus its complex conjugate, and for the spin magnetization, 

- = --(«*F + F*«)-7i(.-C), (20) 

where ( is the equilibrium magnetization of a cluster. 

There exists a large variety of different nanoclusters [1,2,6,35-38], because of which the 
system parameters can take values in a wide range. Concrete examples will be discussed in 
the concluding section. 



3 Approximate Analysis of Spin Dynamics 

The equations of motion (19) and (20) are derived for arbitrary cluster parameters. Their 
solution can be done numerically. But before we pass to numerical investigation, it is useful 
to give an approximate qualitative analysis allowing for the better understanding of physics 
involved. For this purpose, wc shall assume that some of the parameters are small as com- 
pared to the Zeeman and resonator frequencies. 

First, it is possible to show [32-34] that the feedback equation (4) can be represented as 
the integral equation 



H = -Attt] I G{t - t')m^{t') dt' , 
Jo 



(21) 



in which the transfer function is 



G{t) = cos{ujt) — — shi{ujt) 

L U! 



e 



with the shifted frequency 
and the source is 



The resonator feedback field efficiently acts on the cluster only when the resonator natural 
frequency is tuned close to the Zeeman frequency, so that the detuning be small, satisfying 
the resonance condition 

A 

— < 1 {A = uj-ujo) . (22) 

The effective rotation frequency (17) should also be close to which requires that the 
anisotropy frequency (16) be also sufficiently small, 

^ « 1. (23) 

OJ 

We assume that all attenuation parameters are smaller than ou, such that 

^ <Cl , — < 1 . (24) 

UJ CO 

And let us introduce the parameter playing the role of the feedback rate 

which characterizes the attenuation caused by the coupling between the cluster and res- 
onator. 

Under these conditions, the solution to the integral equation (21) can be found by an 
iterative procedure [32-34], which here gives 

/IqH — ih {au — a*u*) , (26) 

where the coupling function, for small detuning, such that | A| < 7, reads as 

77o^ 



a 



72 + A2 



(1 - e-^*) . (27) 



Substituting form (26) into Eqs. (19) and (20) results in the equations 

— = —iilu + asu + u!is — s (au)* , 



dt 



— = 2asw + u!i {u* + u) s — 2sRe (au^) , 



— — —aw — — {u* + u) — ^i{s — Q + Re (au^) , 



(28) 




(29) 



According to inequalities (22) to (24), the variables w and s can be treated as slowly 
varying in time, while w as a fast variable. The existence of two time scales allows us to 
invoke the scale separation approach [30,31,39,40] that is a variant of the averaging techniques 
[41]. Then we solve the equation for keeping there w and s as quasi-integrals of motion, 
which gives 



where Uq = m(0). This solution is to be substituted into the equations for the slow variables 
w and s, with averaging their right-hand sides over time, thus, obtaining the equations for 
the guiding centers of w and s. 

The scale separation approach is a very powerful method for dealing with a system of 
many interacting elements, such as the assemblies of spins [30-34] or of quantum dots [42]. 
However, it requires the use of limitations on the range of the system parameters, such as 
inequahties (22) to (24). While, in our case, we have derived Eqs. (19) and (20) that are valid 
for arbitrary parameter values. The latter equations are not difficult to solve numerically. 
But the present approximate consideration is useful for understanding the following physical 
points. 

The source of creating the resonator feedback field (21) is the motion of the cluster spin, 
which generates the magnetic field (26) acting back on the spin. The back action is absent 
at the initial time. The coupling of the cluster spin and the resonator increases with time 
according to the behavior of the coupling function (27). The transverse spin component (30) 
oscillates with time with the rotation frequency (17). The initial push to the spin oscillation 
is done by the transverse field Bi entering through the transverse frequency (7). This also 
requires that the initial magnetization (12) be nonzero. Generally, there is one more source 
triggering the initial spin motion. These are quantum spin fluctuations that start the spin 
motion even in the absence of a transverse field, as has been shown for multi-spin systems 
[30-33]. Such quantum spin fluctuations are especially important for particles with low spins. 
The input of quantum fluctuations for large spins 5" ^ 1 diminishes as 1/5", in agreement 
with Eq. (13). Therefore, in the case of a single nanocluster, these fluctuations can be 
neglected as soon as the transverse fleld Bi is present. The transverse fleld, that can be 
easily regulated, serves as a convenient triggering mechanism for initiating spin motion. 

4 Numerical Solution of Evolution Equations 

Now we go back to the general evolution equations (19) and (20). For numerical investigation, 
it is appropriate to pass to dimensionless quantities. To this end, we deflne the dimensionless 




(30) 



feedback field 

h^-^ = ^-^. (31) 

«7o 7o 

The effective forces (8) and (18), respectively, become 

/ = i7o/i + , F = i7o/i + - Cc;2W* + ■ (32) 
Instead of the complex variable (11), let us use the real components 

{S^) , , (33) 

when ^ 

?i = - {S-) =x-iy . 

And in the following, let us measure time in units of l/7o- 

With these notations, Eqs. (19) and (20) yield the equations for the spin components: 

— = -ojoy + uoys + uis , (34) 
dy 

— — ujqX — (ub + uj2)xs — hs , (35) 

(JjV 

ds 

— ^ hy - ujix + uj2xy - -fi{s - C) . (36) 
The feedback equation (4) can be rewritten in the form 

d'^h ^ dh 2 7 A /o^,\ 
_ + 2,-+.^/. = 4-. (37) 

We assume that at the starting moment of time the cluster magnetization is polarized along 
the axis z, which implies the initial spin components 

xo = x(0) = 0, ?/o = y(0) = 0, so = s(0) = l. (38) 

And for the initial feedback field we set 

ho = h{0) = , ho = h(0) = , (39) 

where the overdot signifies time derivative. 

The overall time evolution is described by Eqs. (34) to (37), with the initial conditions 
(38) and (39). These equations possess the stable stationary solution 

X* = , y* = ^ , s* = C, h*^h*^0, (40) 

Uo - UJdC 

that is reached in the relaxation time Ti — I/71. The attenuation 71 can be defined by the 
Arrhenius law 

Ea 



71 = 7a exp 



where Ea = hwAS is the anisotropy barrier. At low temperatures, below the blocking 
temperature, the attenuation 71 is exponentially suppressed, so that the relaxation time Ti 
is very long. However the magnetization reversal can be ultrafast due to the action of the 
resonator feedback field. In our calculations, we set 71 = 10~^ (in units of 70) and take the 
initial conditions (38) and (39). 

Figure 1 shows the solutions to the evolution equations (34) to (37), with the initial 
conditions (38) and (39), as functions of time (in units of I/70) for typical parameters 
expressed in units of 70. The resonator feedback field realizes the magnetization reversal 
that can be many orders shorter than the relaxation time Ti. Comparing the behavior in 
the presence of the resonator and in its absence, we see that in the time scale of I/70 the 
influence of 71 is not noticeable at all. 

Figure 2 demonstrates the infiucnce of the resonator attenuation on the magnetization 
reversal. Too short 7 means a large ringing time I/7, when the resonator several times 
exchanges the energy with the cluster, which leads to the oscillations of the magnetization. 
The optimal value of the resonator attenuation is 7 = 1, when the resonator attenuation 7 
coincides with the attenuation 70 characterizing the coupling between the resonator and the 
cluster. 

Figure 3 illustrates the role of the Zeeman frequency. The larger the latter, the shorter 
the reversal time. However too large (Jq leads to many oscillations of magnetization, which 
is not good, if one aims at the stable reversal. 

Figure 4 shows the influence of the transverse field Bi entering through the transverse 
frequency Ui. The larger the latter, the shorter the reversal time. But too large Ui results 
in multiple oscillations, which may be inconvenient for practical purposes. 

Figure 5 demonstrates the role of magnetic anisotropy. The anisotropy frequencies that 
are smaller than the Zeeman frequency do not strongly influence the reversal. But if the 
anisotropy frequencies are larger than the Zeeman frequency, then the reversal is blocked. 

Figure 6 shows the importance of the resonance condition. Large detuning from resonance 
makes the magnetization reversal slower. The larger the detuning, the slower the reversal. 
The reversal is the fastest when the resonance condition cjq = <^ is valid. 

These figures demonstrate that the coupling to a resonant electric circuit results in the 
ultrafast magnetization reversal of a single nanocluster, as compared to its natural relaxation 
time caused by thermal fiuctuations and phonon-assisted tunneling. Estimates for typical 
nanoclusters will be given in the concluding section. 

5 System of Nanoclusters with Dipolar Interactions 

In the previous sections, we have treated the magnetization reversal of a single nanocluster. 
An important question is whether such an ultrafast magnetization reversal could be achieved 
for an ensemble of nanoclusters. The basic difference of the latter case from that of a 
single cluster is the existence of strong dipolar interactions between the nanoclusters. These 
dipolar interactions completely suppress coherent spin motion in dephasing time T2, so that 
collective spin rotation, without a resonant feedback, becomes impossible [32-34,43]. One 
should not confuse the real dipolar interactions between spins with the effective interactions 
through photon exchange of resonant atoms radiating at optical frequencies. The dipolar 
spin interactions dephase spin motion, while the atomic interactions through the photon 
exchange, vice versa, collectivize atomic radiation [44]. Spin motion, of course, also produces 



electromagnetic radiation that, however, is extremely weak and can never collectivize spins in 
time shorter than the dephasing time [33,34,43,44]. Self-organized coherent atomic radiation, 
called superradiancc, is the Dickc effect [45] . The principally different collectivization of spin 
motion by means of a resonator feedback field is what is termed the Purcell effect [26]. The 
Dicke effect for spin systems is impossible, so that the self-organized coherent spin motion 
is admissible only through the Purcell effect [32-34,43,44], which necessarily requires the 
coupling of the spin system with a resonator. 

Let us consider a system of N nanoclusters in volume V, with the density p = N/V . The 
system Hamiltonian reads as 

N N 

H = Y.H^ + 2^H,,. (41) 

Here the first term is the sum of the single-cluster Hamiltonians 

4- = -/XoB • S,- - D{S'^f + L>2(5J)' , (42) 

with the cluster spin operators Sj and the index j = 1,2, ... ,N enumerating the nan- 
oclusters. Aiming at studying the role of the dipolar interactions, we keep in mind similar 
nanoclusters with close anisotropy parameters and spins S. The dipolar interactions are 
characterized by the Hamiltonian parts 

H^,=Y.^ff^^^f ^ (43) 
a/3 



with the dipolar tensor 



jja/3 _ /fo 
ij 



- 3<<.) , (44) 



in which 



II ^ij 



The resonator feedback field is described by the same Eq. (4), but with 



TV 



To write the equations of motion in a compact form, we introduce several notations. The 
dipolar terms are combined into the variables 



(46) 



where 



a,,^D-, h,^\{Ul^-iy^-'MD^) , c,,^l{Dtf-tDy;) . (47) 



In the evolution equations, as is well known, there arise pair spin correlators that need to 
be decoupled for obtaining a closed system of differential equations. For different clusters, 
we use the semiclassical decoupling 

{srsf)^{sr){sf) (i^j), (48) 

supplemented by the account of quantum spin correlations yielding the appearance of the 
dephasing term 72. And the spin correlators for different components of the same cluster, 
similar to Eq. (13), are decoupled as 

{SfS^ + S^S^) = (2 - {S-) {S^) (a 7^ /3) , (49) 

in order to retain the correct limiting expressions for spin one-half and large spins [32-34]. 
The angle brackets, as earlier, imply statistical averaging over the initial statistical operator. 
This type of spin decoupling will lead to the appearance of the expressions 

N N 

As we see, the consideration of an ensemble of magnetic nanoclusters with dipolar in- 
teractions is essentially more complicated than that of a single nanocluster, treated in the 
previous sections. 

6 Magnetization Reversal in Ensemble of Nanoclusters 

The quantities of interest are the average spin components, for which it is convenient to 
introduce, instead of Eqs. (11) and (12), the relative transverse component 



and the relative longitudinal component 

1 ^ 

i=i 

The effective force acting on a cluster spin, instead of Eq. (18), now is 

F = - ^ l^oH + UJ1 + I UJ2U* + C ■ (53) 

Writing down the equations of motion for the spin operators and averaging them [46], we 
come to the evolution equations for the mean spin components (51) and (52) in the form 

nil n ^ 1 

^ = -^(^] + eo - iV2)u + Fs, ^ = - 2 («*^ + P*^) - 7i(^ - C) , (54) 

where 

r2 = 72(l-«^), ^.^^^ph^lS. (55) 



Equations (54) for a nanocluster system replace Eqs. (19) and (20) for a single cluster. 
These equations are to be considered together with Eq. (4) for the resonator feedback field, 
with ^ 

= 2 Pl^oS{u* + u) . (56) 
Instead of the feedback rate (25) for a single cluster, for N clusters, we now have 

7o(iV) = ttTV = Trr/p^l^ • (57) 

And instead of the couphng function (27), we get 

ot{N) = gi2{l - As) (1 - e-^*) , (58) 

with the dimensionless coupling parameter 

_ 7a;o7o(A^) /.on 
« = ^;h?TA5) (59) 

and the effective anisotropy parameter 

^^a^^2a^^+a^_ (60) 
ojo 2a;o 

Notice that, for nrj ~ 1, the feedback rate (57) is of order of 72. Under good resonance, with 
a small detuning A ~ 0, the coupling parameter (59) is g ^ ujq/j. 

The system of equations (54) and (21) can be solved numerically, either directly, as has 
been done for magnetic molecules [47], or invoking the averaging techniques [39,40], when fast 
oscillations are averaged out, so that the resulting curves are smoothed. Both these methods 
give close results. The averaging techniques provide more physically transparent description 
of the initial stage of spin motion, showing that this motion starts with stochastic spin 
fluctuations caused by nonsecular terms of dipolar interactions. Thus, dipolar interactions 
play at the initial stage the positive role of a triggering mechanism initiating spin motion, 
while at the later stage they play the negative role, by destroying the coherence of spin 
rotation. 

In Fig. 7, we show the results of numerical calculations, involving the averaging tech- 
niques, for the time dependence of the reduced spin polarization (52) for different system 
parameters. Time is measured in units of 7^^ and all frequencies, in units of 72. The realistic 
value for the anisotropy parameter (60) is taken as A = 0.1. It is worth emphasizing that 
for 72 > cji, dipolar interactions initiate spin dynamics so that coi plays a minor role. In 
order to stress the triggering role of the dipolar interactions, the transverse frequency is set 
to zero, ui — 0. Analogously to Eqs. (38), the initial conditions are 

Wo = w{0) = , So = s(0) = 1 , 

where w = |Mp = a;^ + y^. The figure shows that ultrafast magnetization reversal happens 
also in a system of nanoclusters interacting through dipolar forces. Even more interesting is 
the fact that the dipolar interactions play the role of a triggering mechanism starting spin 
dynamics. The magnetization reversal is realized during the time of order l/g^2 = T2/g, 
which, for g > 1 is shorter that the dephasing time T2 = 1/72- Hence it is feasible to find the 
system parameters, when dipolar interactions do not disturb the coherence of spin motion, 
provided the sample is coupled to a resonator. Without the latter, the Purcell effect does 
not exist and the ultrafast magnetization reversal is impossible. 



7 Discussion 



We have suggested a method for reahzing an ultrafast magnetization reversal of nanoclusters. 
The possibihty of such a fast reversal is important for a number of applications, e.g., for 
the functioning of various magneto-electronic devices, spintronics, magnetic recording and 
storage, and other information processing techniques [48-50] . The idea of the method is based 
on the coupling of the nanocluster with a resonant electric circuit. This is easily achievable by 
placing the nanoclusters inside a magnetic coil. Then the motion of the nanocluster magnetic 
moment produces electric current in the circuit, which creates magnetic field acting back on 
the cluster magnetization. This feedback field of the resonator accelerates the magnetization 
reversal. The reversal time can be made many orders shorter than the natural relaxation 
time. 

First, we have considered a single nanocluster, which makes it possible to avoid compli- 
cations due to distributions of particle sizes, shapes, spin values, and so on, which could arise 
in the case of an assembly of many nanoclusters. In the latter case, the basic complication 
is the necessity of taking into account dipole interactions between the clusters. All these 
additional problems are avoided when dealing with a single cluster. 

The case of an ensemble of nanoclusters, interacting through dipolar forces is also an- 
alyzed. The ultrafast magnetization reversal is feasible for this case as well. The reversal 
occurs during the time shorter than the dipole dephasing time, because of which dipolar 
interactions do not destroy coherent spin motion that is responsible for the ultrafast rever- 
sal. Even more, dipolar interactions are useful at the initial stage, when they trigger spin 
dynamics. 

To give the feeling of typical values for the characteristic parameters, let us make esti- 
mates for some nanoclusters. Actually, the family of magnetic nanoclusters is very wide and 
these can display rather different properties [1-3,51-54]. To be concrete, let us keep in mind 
the values typical of Co, Fe, and Ni nanoclusters. The coherence radius for these clusters, 
below which they are in a single-domain state and can display coherent rotation of magne- 
tization is Rcoh ~ 10 nm. The standardly formed clusters have the radii i? ~ 1 — 3 nm. The 
corresponding cluster volume is Vi ~ lO"^'^ cm^. A cluster contains about Ni ~ 10^ atoms. 
The atomic density in a cluster is pi ~ 10^'^ cm~^. The cluster spin is proportional to the 
number of atoms in the cluster, hence S* ~ 10^, that is, the magnetic moment is of the order 
10^ Hb- The magnetic anisotropy parameters (2) are Ki ~ K2 ~ 10^ erg/cm^. The fourth- 
order anisotropy is much smaller, ^^^4 ~ 10^ erg/cm^. This, for the anisotropy parameters of 
Hamiltonian (1), gives D D2 ^ 10"^*^ erg. And for the fourth-order anisotropy, this would 
make ~ 10~^^ erg. 

These values, for the anisotropy frequencies (14) to (16) yield uo ~ ~ ~ 10^*^ 
Hz. This corresponds to the anisotropy field Ba = oja/'^s ~ 10^ G. The Zeeman frequency, 
for the magnetic field ~ 1 T is a;o ~ 10^^ Hz. Note that the present day facilities allow 
for the generation of magnetic fields as high as about 100 T [55] . The feedback rate (25) is 

70 ~ 10^*^ s~^. This rate provides the reversal time trev ~ l/o'o ~ 10""*"° s. 

The blocking temperature, below which thermally activated reversals are exponentially 
suppressed is Tb ~ 10 — 40 K. The typical prefactor in the Arrhenius law is 7^ ~ 10^ — 10^^ 
s"-*^. The anisotropy energy barrier in the Arrhenius law is Ea ~ lO"-*^^ erg. This gives the 
anisotropy temperature EA/ks ~ 100 K. The resulting relaxation time Ti = I/71, below 
the blocking temperature, is rather long. Thus, even at the temperature T = 10 K, we have 

71 ~ 10~^ s. At the temperature T = 5 K, one has Ti ~ 0.1 s. And for temperature T — 1 



K, the thermal reversal time is astronomically large, being Ti ~ 10 s. But, coupling the 
nanocluster to a resonator, produces a very short reversal time tj-ev ~ 10~^° s, independently 
of the value Ti. The reversal time could be made shorter by choosing the appropriate types 
of nanoclusters and resonator properties. 

As has been explained above, the spin relaxation in the system of nanoclusters, which 
would be caused by the photon exchange between different spins is negligible [33,34,43,44]. 
The corresponding radiation width is 

For the typical density of nanoclusters p ~ 10^° cm~^, the spin of a cluster S ~ 10^, and 
frequency u) ~ 10^^ Hz, we get ^rad ~ 10"^ s~^, which is much smaller than the dipolar 
width 72 ~ 10^° s~^. Therefore the relaxation time, due to the photon exchange between 
spins, trad = '^/irad ~ 10® s ~ 10 ycars, is enormously larger than the dipolar dephasing 
time T2 ~ 10~^° s. This confirms that the photon exchange mechanism plays no role in the 
spin relaxation, that is, the Dicke effect for spin systems does not exist. But the ultrafast 
magnetization reversal is completely due to the Purcell effect. 

The typical density of an ensemble of clusters is p ~ 10^° cm~^. With the natural dipolar 
width 72 ~ 10^" s~^, the dephasing time is T2 ~ 10~^° s. As is seen in Fig. 7, the reversal 
time can be an order shorter than the dephasing time, being t^ev ~ 10^^^ s. 

Concluding, by coupling nanoclusters to a resonant electric circuit, it is possible to realize 
ultrafast magnetization reversal for single nanoclusters as well as for assemblies of nanoclus- 
ters. For such nanoclusters as Fe, Co, or Ni, the reversal time can be made as short as 10~^^ 
s. 
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Figure Captions 



Fig. 1. Solutions to the evolution equations: (a) x = x(t); (b) y — y{t); (c) s = s{t); 
(d) h — h{t) for the parameters ouq — uj — 10, ujd — oui = U2 — 1, ^ — — 10"^. The 

parameters are measured in units of 70 and time, in units of I/70. To emphasize the role of 
the resonator feedback field, the solutions in the presence of the resonator (solid lines) are 
compared with those for the case of no resonator (dashed lines). 

Fig. 2. Role of the resonator attenuation. Magnetization as a function of time for the 
parameters uq = u = 10, ujd = uji = uj2 = I, and varying resonator attenuation: (a) 7 = 0.1; 
(b) 7 = 1 (solid line) and 7 = 10 (dashed line). 

Fig. 3. Role of the Zccman frequency. Magnetization as a function of time for the 
parameters u^, = Ui = UJ2 = 1, ^ = 1, and varying Zeeman frequency: (a) a;o = cu = 1 (solid 
line) and cjq = = 10 (dashed line); (b) cuo = ou = 100. 

Fig. 4. Role of the triggering field. Magnetization as a function of time for the parame- 
ters ojq = u = 10, ud = U2 = I, and varying triggering field: (a) cui = 0.001 (solid line) and 
cji = 1 (dashed line); (b) cui = 10. 

Fig. 5. Role of the anisotropy. Magnetization as a function of time for the parameters 
cuq = UJ = 10, cui = 1,7= 1, and varying anisotropy frequencies = a;2 = 1 (dashed-doted 
line); (jJd = a;2 = 10 (solid line); cud = 0^2 = 15 (dashed line). 

Fig. 6. Role of the resonance. Magnetization as a function of time for the parameters 
cuq = 10, cu/j = cui = UJ2 = 1, ■J = 1, and varying detuning from the resonance, with a; = 1 
(sohd line); a; = 10 (dashed-dotted fine); cu — 20 (dashed line). 

Fig. 7. Magnetization reversal in the system of nanoclusters with dipolar interactions. 
The system parameters are 71 = 10~^, A = 0, ^4 = 0.1, cJi = 0. Time is measured in units 
of T2 = 1/72 = 10~^° s, and the frequencies, in units of 72. Other parameters are: 7 = 10, 
uq — u! — 1000, g — 100 (sohd fine); j — 1, loq — u — 100, g — 100 (dashed-dotted line); 
7 = 10. uo = UJ = 100, g = 10 (dashed line). The shown functions of time are: (a) coherence 
intensity w{t); (b) reduced magnetization s{t). 
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Figure 1: Solutions to the evolution equations: (a) x = x(t); (b) y = y{t); (c) s = s(t); 
(d) h = h(t) for the parameters uq = u = 10, ud = i^i = 1^2 = 1, 7 = 1, 7i = 10~^. The 
parameters are measured in units of 70 and time, in units of 1/70- To emphasize the role of 
the resonator feedback field, the solutions in the presence of the resonator (solid lines) are 
compared with those for the case of no resonator (dashed lines). 




Figure 2: Role of the resonator attenuation. Magnetization as a function of time for the 
parameters Uq = u = 10,Ud = Ui = U2 = I and varying resonator attenuation: (a) 7 = 0.1; 
(b) 7 = 1 (sohd hne) and 7 = 10 (dashed hne). 




Figure 3: Role of the Zeeman frequency. Magnetization as a function of time for the param- 
eters Ud = uji = <-02 = I, 'J = I and varying Zeeman frequency: (a) Uq = u = 1 (solid line) 
and uo = u = 10 (dashed line); (b) uo = u = 100. 




Figure 4: Role of the triggering field. Magnetization as a function of time for the parameters 
ujq = uj = 10, ujd = c<;2 = 1, and varying triggering field: (a) Ui = 0.001 (solid line) and 
cji = 1 (dashed line); (b) ui = 10. 




Figure 5: Role of the anisotropy. Magnetization as a function of time for the parameters 
luq = u = 10, Ui = 1, •y = 1, and varying anisotropy frequencies un = U2 = 1 (dashed-doted 
line); Ud = U2 = 10 (sohd line); cod = c<;2 = 15 (dashed line). 
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Figure 6: Role of the resonance. Magnetization as a function of time for the parameters 
ujo = 10, ud = wi = ci;2 = 1, 7 = 1, and varying detuning from the resonance, with a; = 1 
(sohd hne); a; = 10 (dashed-dotted hne); u = 20 (dashed hne). 




Figure 7: Magnetization reversal in the system of nanoclusters with dipolar interactions. 
The system parameters are 71 = 10"'^, A = 0, A = 0.1, Wi = 0. Time is measured in units 
of T2 = 1/72 = 10"^" s, and the frequencies, in units of 72. Other parameters are: 7 = 10, 
ujq = u = 1000, g = 100 (sohd line); •y = 1, Uq = u = 100, g = 100 (dashed-dotted line); 
7 = 10, Wo = u; = 100, g = 10 (dashed line). The shown functions of time are: (a) coherence 
intensity w(t); (b) reduced magnetization s(t). 



